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y-^ . Abstract — The study of sterile neutrinos has recently acquired a different flavor: being now excluded 

' as the dominant solution for the solar or atmospheric conversions, sterile neutrinos, still attractive for 

many other reasons, have thus become even more elusive. The present relevant questions are: which 
subdominant role can they have? Where (and how) can they show up? Cosmology and supernovae 
turn out to be powerful tools to address these issues. With the most general mixing scenarios in mind, 
' I present the analysis of many possible effects on BBN, CMB, LSS, and in SN physics due to sterile 

CSJ i neutrinos. I discuss the computational techniques, present the state-of-the-art bounds, identify the still 

allowed regions and study some of the most promising future probes. I show how the region of the 
LSND sterile neutrino is excluded by the constraints of standard cosmology. 
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O ■ 1 Introduction 



The study of sterile neutrinos (namely: additional light fermionic particles that are neutral under all Standard Model 
gauge forces, but can be a non-negligible ingredient of our world through their mixing with e, /i, r neutrinos) has 
recently acquired a different flavor: the established solar and atmospheric anomalies seem produced by oscillations 
among the three active neutrinos |2| [S] so that their explanation in terms of oscillation into a Vs state (which was 



o 

. viable and fairly popular up to a few years ago) is believed to be now ruled out as the dominant mechanism. 

This means that the relevant questions concerning sterile neutrinos nowadays have become: •which is the 
subdominant m\e still possible for sterile neutrinos in solar and atmospheric neutrinos? •where can we detect the 
' effects of a sterile neutrino? i.e. which are the most sensitive experiments (in astrophysics, cosmology or man- 

made set-ups) in which sterile neutrinos can be discovered? •how can we detect the effects of a sterile neutrino? 
i.e. which are the signatures of its presence? 

To answer these questions, the investigation on sterile neutrinos requires a more extensive and deep approach. 
Indeed, for instance, most of the previous analysis used to consider only the peculiar oscillation pattern that gives 
the simplest physics: the initial active neutrino \va) {ve in the case of solar neutrinos, in the case of the atmo- 
spheric ones) oscillates into an energy-independent mixed neutrino cos9s\v'^) + sm 6 s\vs). This of course leaves 
unexplored the largest part of the parameter space. Sometimes, moreover, previous analysis used to neglect for 
simplicity the mixing among active neutrinos (this is the case of most studies on the sterile effects in cosmology or 
in supernovae). Such a mixing is now established and its parameters are reasonably pinned down. 

These considerations motivate the analysis performed in |fTl and presented here, which considers and includes: 



a. any possible l'e,^l,T — mixing pattern; 

b. the established v,, — v^,t, — Vr mixings; 

c. all possible neutrino sources and contexts (cosmology (BigBang Nucleosynthesis-BBN, Cosmic Microwave 
Background-CMB, Large Scale Structures-LSS), astrophysics (the Sun, supernovae-SNe...), atmospheric 
neutrinos, reactor and accelerator experiments...), 

*)Based on the Proceedings for the lO"^ International Symposium on Particles, Strings and Cosmology (PASCOS '04), 16-22 August 2004, 
Northeastern University, Boston, MA, USA and for the XVI Incontri suUa Fisica delle Alte Energie (IFAE), 14-16 April 2004, Torino, Italy. 



The purpose is to set the state-of-the-art bounds on the active-sterile mixing parameters (in each context separately 
and then in a combined way) and to identify the most promising future probes. Since the different cosmological 
quantities and the observables of SN physics turn out to be very important to this aim, investigating complementary 
regions and allowing the application of techniques that extend to other fields, these are the contexts on which I focus 
in these Proceedings. 

Before going to the details, let's stress that now that the hunger of sterile neutrinos for the solar and atmospheric 
anomalies is over, nevertheless sterile neutrinos are even more attractive for at least two sets of reasons. 
First, from a top-down point of view, sterile fermions that are naturally light or cleverly lightened arise in many 
theories that try to figure out what is beyond the Standard Model. To begin, slightly beyond the context of the 
SM, the right handed neutrino is an obvious candidate, to complete the lepton sector in similarity (symmetry?) 
with the quark one. In this case, actually, three states (one per family) would be natural. More broadly, several 
(SuSy-/GUT-/string-/ED- inspired) SM gauge singlets line up awaiting for consideration (axino, branino, dilatino, 
familino, Goldstino, Majorino, mirror fermion, modulino, radino...) yj. The origin and the load of information of 
any of these particles can be very different, but from an effective point of view we simply need to parameterize 
their mixing with the active neutrinos in terms of mixing angles and Am^ (see below) to include them all in the 
analysis. Independently from the specific model, the discovery of a new light particle would be of fundamental 
importance and deserves to be investigated per se. 

Second, from a more phenomenological bottom-up perspective, sterile neutrinos are repeatedly pointed as a pos- 
sible explanation for several puzzling situations in particle physics, astronomy and cosmology. For instance, they 
have been invoked |5l to account for the origin of the pulsar kicks, to constitute a Dark Matter candidate, to explain 
(via their decay) the diffuse ionization of the Milky Way, to help the r-process nucleosynthesis in the environment 
of exploding stars, to interpret the slightly too low Argon production in the Homestake experiment... and notori- 
ously to explain the LSND claimed evidence of oscillations |6l. Any one of these puzzles, in general, points to 
specific sterile neutrinos (i.e. with a specific mixing pattern with the active ones) and it is therefore worthwhile to 
explore them in an extensive way. 

2 Four neutrino mix 

In absence of sterile neutrinos, we denote by U the usual 3x3 mixing matrix that relates neutrino flavor eigenstates 
i^e,fj..T to active neutrino mass eigenstates 1^1,2.3 as i^i = UuVi- The extra sterile neutrino can then mix with 
a mixing angle 6s with an arbitrary combination of active neutrinos, which we identify by a complex unit 
3 -vector n 

n-v^ neVe + n^v^ + UrVr = "i^^i + «2i^2 + n^v^ {ui = Uuni). (2) 

So, in particular, the 4"^ mass eigenstate is given by 1/4 = i^g cos 6s + nii/( sin 6s. We allow all possible values of 
its mass ('the mass of the sterile neutrino', in the small mixing case). 

Such a formalism is completely general and covers of course all the possible mixing patterns. We need to choose, 
however, some intuitive limiting cases to present the results in the following: 

• Mixing with a flavor eigenstate (fig.^): the sterile neutrino oscillates into one of the active flavors {n-v = 
vi with ^ = e or ^ or t). Therefore there are 3 different active-sterile Arri^ (which cannot all be smaller 
than the observed splittings Am^m^ g^^j^^, see figure). 

• Mixing with a mass eigenstate (fig.E')- the sterile neutrino oscillates into a matter eigenstate, that consists 
therefore of mixed flavor {fi-v ^ Vi, with i = 1 or 2 or 3). There is one single Am^, which can be arbitrarily 
small. 

^ ) In this way n and ds , together with the four neutrino masses, simply reorganize in a more intuitive way all the parameters of the most genetic 
4x4 Majorana neutrino mass matrix, which contains 4 masses, 6 mixing angles and 6 CP-violating phases, of which 3 affect oscillations. The 
4x4 neutrino mixing matrix V that relates flavor to mass eigenstates as i^e,n,T,s = V ■ 1^1,2,3,4 is expressed in this parameterization by 

, , _ / 1 — (1 — cosSs)?!* ® n sin6*sn* \ ( U Q\ 

\ -sinSsn cos6»s ) ^ \ f) l) ^ ' 

Commonly, V is instead parameterized asV = R-nR2iRiA ■ U23Ui:jUi2 when studying sterile mixing with a flavor eigenstate. Here Rij 
represents a rotation in the ij plane by angle 6ij and Uij a complex rotation in the ij plane. 614 or Ue4 give rise to i/g /us mixing, 824 or U/^i 
to i^fi/us mixing, and 934 or Uri to Vr jva mixing. When studying the mixing with a mass eigenstate, the parameterization of V is changed to 
y = U23U13U12 ■ -R34 -R24 ^14 ■ Now 9n gives rise to Vijv^ mixing and so on. Our parameterization is more convenient because n already 
encrypts in a natural way the information on which active states mix with the Vs, while 0s simply expresses the size of the mixing. 
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Figure 1: Basic kinds of four neutrino mass spectra. Left: sterile mixing with a flavor eigenstate (v^ in the picture). 
Right: sterile mixing with a mass eigenstate ( 1^2 in the picture). 

In the case of mixing with a matter eigenstate, we will consider also the situation in which the (mostly) sterile state 
is lighter than the (mostly) active state with which it mixes (imagine the 1/4 state lowered below 1^2 in figure^): it 
is represented by the portion of 9s > it /A in our plots. 

We assume that active neutrinos have normal hierarchy, Am23 > 0. Finally, we assume 9is = 0. We verified that 
using 6*13 ~ 0.2, the maximal value allowed by present experiments, leads to minor (in some cases) or no (in other 
cases) modifications, that we do not discuss. 

3 Sterile Neutrinos in cosmological sauce 

Generalities: The Early Universe can be a powerful laboratory for neutrino physics, and therefore in particular for 
the physics of sterile neutrinos. The fundamental reasons for this basic fact are simply listed: 

• (light) neutrinos are very abundant (namely "as abundant as photons") for a long period of the evolution of 
the Universe, keeping thermal equilibrium until T ~ few MeV; 

• in a Friedman-Robertson- Walker standard cosmology, the total energy density is a crucial parameter that 
sets the expansion rate of the Universe; since that energy is predominantly stored in the relativistic species, 
namely electrons, positrons, photons and all species of neutrinos, for T ~ 100 MeV 1 MeV, it is evident 
that the relative abundance of neutrinos (e.g. increased by the presence of additional states) is a very relevant 
quantity; 

• the early plasma is so dense that neutrinos are initially trapped and undergo peculiar matter effects while the 
density decreases as a consequence of the expansion; 

• the detailed balance of the different species of neutrinos among themselves can also be important for pro- 
cesses that distinguish flavor: for instance, the ly^ density affects the n ^ p conversion and therefore is 
imprinted in the primordial ratio of n/p that we read today (see below). 

We have access to several different windows during the history of the Universe. From them we get a sensitivity to 
neutrino properties (masses, oscillation parameters...) that is nowadays competitive with direct measurements and 
even offers a brighter prospective of improvements in the near future, making the study of quantitative neutrino 
cosmology worthwhile. 

BBN: BigBang Nucleosynthesis occurs at T ~ (1 0.1) MeV and describes the era when the light elements 
were synthesized |[7j. Given a few input parameters (the effective number iV^ of thermalized relativistic species, 
the baryon asymmetry ub /n-y ^ rj, and possibly the i^t j vi lepton asymmetries) BBN successfully predicts the 
abundances of several light nuclei. Today r\ is best determined within minimal cosmology by CMB data to be 77 = 
(6.15 ± 0.25)10~^° jSj. Thus, neglecting the lepton asymmetries (see below), basically one uses the observations 
of primordial abundances to test if = 3 as predicted by the SM-^'. This is what is often done 1101 . using state- 
of-the-art, publicly available codes fl Ij. We need however to do something slightly more refined: indeed, Ny 
is an effective parameter that sums up and does not distinguish the relative contributions of the various neutrinos 

^'For the sake of precision, relaxing the hypothesis of instantaneous neutrino freeze-out and carefully including the partial neutrino reheating 
from e+e~ annihilations, the related spectral distortions and tinite temperature QED small effects, the SM prediction is actually N^, ~ 3.04 1^. 
The deviations due to any exotic phenomenon (including sterile neutrinos) go on top of this. 
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(3 active and 1 (or more) sterile), which can instead be different according to the mixing patterns; moreover, in 
general, the neutrino populations can have a peculiar behavior in time (temperature) which would be concealed by 
the use of N^,. 

For these reasons, in our computation we use as variables the four neutrino densities '^^ p^^ , pi, , p^^ , p^_^ : we 
follow their evolution with the temperature during the whole period of BBN, for each possible choice of the sterile 
mixing parameters. Starting (conservatively) from a zero initial abundance of the sterile neutrino, at a certain 
point oscillations start producing it |,12 |. This essentially happens when the plasma effects/thermal masses for the 
(active) neutrinos cease to be dominant compared to the vacuum masses, as the Universe expands and cools; at 
what point precisely (and how efficiently) the production occurs is something which is determined by the sterile 
mixing parameters. In the meanwhile, other cosmological processes occur, as standard: at T few MeV neutrinos 
freeze-out (i.e. loose thermal equilibrium with the bath, but they still take part in the n ^ p reactions, see below); 
at T ~ 1 MeV e+e~ annihilate etc... 

More precisely, we follow the time evolution of the full 4x4 density matrices g and q (written in the flavor basis), of which the four neutrino 
densities above constitute the diagonal. In absence of neutrino asymmetiies, the neutrino and antineutrino sectors decouple and proceed 
identically, so we focus on the neutrinos for definiteness. The kinetic equations for such matrices must take into account (i) the vacuum 
oscillations (active-active and active-sterile), (ii) the matter effects in the primordial plasma, (iii) the i/e ue scattering reactions and the 
uu « ee annihilation reactions. They read I12II13|IT^ 

Hm is the Hamiltonian in matter, composed by the vacuum Hamiltonian in the flavor basis and the matter potentials Vi for each flavor: these 
consist of the thermal masses for the neutrinos in the primordial plasma, rapidly decreasing with T. The usual MSW potential is in this case 
subdominant because the plasma is charge symmetric. 

Hm = \vdia.g{mi,m'^,ml,ml)V^ + E^dag{Ve,V^,VT,0)\ (4) 

= -^^^^G^i^C^TUose^e,,...) (5) 

Vs = 

Concerning the reaction part, it can be shown that one has to use: in the equations for the off-diagonal components of g, Ftot ~ 3.6 G% 
for Ue and Ftot ~ 2.5 Gp for i^^j.t; in the equations for the diagonal components, Fann ~ 0.5 Gp for Ue and Fann ~ 0.3 Gp 
for i^(j,T. q"'^ = diag(l, 1, 1, 0) is the equihbrium value of the density matrix to which the reactions tend. The neutrino freeze-out occurs 
when these F's are overwhelmed by the expansion of the Universe. 

The determination of ^ is in principle quite involved, since we need to keep track of the several phenomena that go on in the range T ~ 1 MeV 
which is under examination. In particular, we want to include the possible extra degrees of freedom (the sterile neutrinos) that are produced by 
the oscillations and we do not want to neglect the heating due to e+e^ annihilations. In first approximation, we can use the standard expression 
T = — (piyjot ) where H contains the (temperature dependent) total energy density, including in particular that of all neutrinos. 

After solving the neutrino densities evolution with temperature, we can study the relative neutrons/protons abundance, which is the all-important 
quantity for the outcome of primordial nucleosynthesis, since those are the building blocks of the nuclei that are going to be formed and 
essentially all neutrons are incorporated into some light element in the process: the neutron abundance at the moment that the synthesis begins 
practically fixes the proportions of all the products, n/p evolves according to 

dT dr Tin 
at al rin + rip 

Tp^n is the total rate for all the p n reactions (n — ♦ p e~ Ve,n — ♦ p , n e+ — > p Pg) and F„^p for the inverse processes. They 
depend on the p^^ and p/^^ densities computed above. At this point, it is apparent how the production of sterile neutrinos can enter the game 
by 

A. entering in Putat ^nd thus increasing the Hubble parameter H i.e. the expansion rate; 

B. modifying the Fp^„, r„_p rates directly, if the Ve, Ve population is depleted by oscillations. 

With the value of n jp in hand, finally a network of Boltzmann equations describes how electroweak, strong and electromagnetic processes 
control the evolution of the various nuclei: p, n, D, T, ^^He, *He,. . . 

For the purposes of the CMB and LSS bounds discussed below, the neutrino densities at the time of recombination and today are needed. These 
are simply given by the final outputs (i.e. for T <C 0.1 MeV) of the kinetic equations. 

We assume a vanishing or negligibly small neutrino asymmetry Allowing an unnatural large (namely, 
much larger than the corresponding asymmetry 77 in baryons), besides affecting non trivially the evolution (15|, 
means adding an extra relevant parameter, so that one could essentially conceal any sterile effect, at least as long 
as observations will not be able to break the degeneracy. 



^^The densities are intended as relative to the photon one, so that p,yj G (0, 1). 
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Figure 2: Some recent experimental results in the determination of the primordial ^He abundance Yp. The bars cor- 
respond to the claimed Icr errors. The PDG recommended value adds an estimate of the systematic uncertainties 
( ±0.005j. References are in M6V . 

At the end of the process, we can compute how the hght elements abundances are modified and compare them to 
the observations. We focus on the ^He abundance, which is today the most sensitive probe of sterile effects, but 
we also study the Deuterium abundance, which has brighter prospects of future improvements. The observational 
determinations of both quantities are plagued by controversial systematic uncertainties. For instance, Fig|2lcollects 
some of the most recent results for "^He flSl . A similar situation, with perhaps less controversy and more overall 
uncertainty, applies to the case of Deuterium 1171 . Conservative estimates are 



Yp = 0.24 ±0.01, 



Y^/Yn = (2.8 ±0.5) 10-5, 
where Yx = nx / ns and Yp is the traditional notation for F4He- 



(7) 



At this point, for ease of presentation, we can even convert the values of the computed primordial abundances back 
into effective numbers of neutrinos, N^^'^ and A^^. For arbitrary values around the SM value of 3, BBN codes 
predict fTl llOl the following relations 



Yp ~ 0.248 + 0.0096 In 
Yb/Yh ~ (2.75 ±0.13) 10" 
The observed abundances of eq.0 then translate into 



+ 0.013(Ar;"^ -3), 



6.15 10-1" 
5 1 + 0.11 (jyp - 3) 
(7]/6.15 10-i")i-6 



N^^'' ~ 2.4 ±0.7, 
iVP ~ 3 ± 2. 



(8 a) 
(8b) 

(9) 



LSS: Neutrinos can also be studied looking at the distribution of galaxies. The connection lies at the time of the 
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formation of the anisotropics in the primordial plasma that were the seeds for the formation of the Large Scale 
Structures (which took place much time later). The point is that the neutrinos, relativistically traveling through the 
plasma (from which they were decoupled) until their mass was of the order of the temperature, had the effect of 
smoothing the anisotropics, i.e. they caused a suppression in the power spectrum of the galaxies that is measured 
today. Qualitatively, light neutrinos traveled relativistically for a long period and therefore delayed the formation 
of structures characterized by a scale smaller than that of the horizon at the time they became non-relativistic. The 
more massive the neutrinos are, the earlier they became non-relativistic, the smaller the scale of the horizon was at 
that time, inside which the perturbations were smoothed, the more suppressed are the large momenta of the LSS 
power spectrum. 

In formulae, the effect is usually expressed in terms of the quantity which is related to the sum of the neutrino 

masses 

93.5 eV 

where m is the 4x4 neutrino mass matrix and g is the 4x4 neutrino density matrix (at late cosmological 
times), discussed and computed above. In a standard case, the numerator corresponds to J2 '^i^i ■ I" general, the 
determination of fi^ depends on priors and on normalisations, possibly fixed by the CMB spectrum. As a rule of 
thumb, the present bound 1181 and the future expected sensitivity 1191 can be summed in 

present: QJi^ <10-^ (11) 

future: n^h^<10~^. (12) 



CMB: Finally, neutrinos can be studied through the pattern of the CMB anisotropics measured by WMAP (and 
other experiments). They affect the CMB anisotropics in various ways |20]; the all important quantity is their 
contribution to the relativistic energy density 

Tr[e] = C prci (13) 

(where again g is the 4x4 neutrino density matrix (at late cosmological times), discussed and computed above), 
straightforwardly parameterized in terms of an effective number of neutrinos 



Prol — P-f 

Global fits at the moment imply IZTl 
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4/3 



^CMB 



(14) 



iV™^«3±2 (15) 



somewhat depending on which priors and on which data are included in the fit. Future data might allow a better 
discrimination. 



3.1 Results 

The results are collected in Fig. |3l We plot the effective numbers Ni, of neutrinos which translate the physical 
observables (the ^He and D abundances and the energy density in neutrinos at recombination), as dictated by 
eqs. ^ and eq. and the value of the quantity il^h^, defined by eq. ilQ\ . We shade the regions that correspond 
to Yp > 0.26 (i.e. A^^^^ > 3.8) or il^h^ > 10^^ and are therefore 'strongly disfavoured' or 'excluded' (depending 
on how conservatively one estimates systematic uncertainties) within minimal cosmology. The other lines indicate 
the sensitivity that future experiments might reach. 

In order to qualitatively understand these precisely computed results, it is useful to begin with the case of mixing 
with a matter eigenstate 1^1,2,3 (upper row of Fig.|3ji, and to consider first the line which corresponds to N^^^ 
(blue dashed line): indeed, this quantity is sensitive only to the total number of neutrinos, and not to the specific 
flavor (j^e, or i^r) which mixes with i/g- In the region above the lines, the production of sterile neutrinos via 
oscillations is efficient and populates to some extent the sterile species. The slope of the lines can be reproduced 
by simple analytical estimates 1121 . and is also intuitive: loosely speaking, at small Am^ the production starts late 

Actually, the naive A^^ discussed at the beginning of the Section is usually defined as equal to this N^^^ . It should be clear from the 
above discussion that, instead, N^^" and stand for two other different quantities. The three of them come to coincide only in the limiting 
case in which the only effect of the extra degrees of freedom is a contribution in the total energy density. 
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Figure 3: Cosmological effects of sterile neutrino oscillations. We collect four different signals. The continuous 
red line refers to the '^He abundance (we shaded as 'strongly disfavoured' the regions where its value corresponds 
to NJ^° > 3.8 i.e. Yp > 0.26j, the purple dotted line to the deuterium abundance, and the dashed blue line to 
the effective number of neutrinos at recombination. We plotted isolines of these three signals corresponding to an 
effective number of neutrinos Ni, = 3.2 and 3.8. The precise meaning of the parameter in the three cases is 
explained in the text. The upper (lower) dot-dashed orange lines corresponds to il^li^ = 10^^ (10^'^)," we shaded 
as 'strongly disfavoured' by the data the regions where Vt^h? > 10~^. 

and needs a large mixing angle to be efficient enough; at large Ajti^, on the contrary, even a small mixing angle 
(i.e. a small rate of production) gives rise to a significant amount. Effects are larger at Os > 7r/4 (i.e. tan^^s > 1) 
because this corresponds to having a mostly sterile state lighter than the mostly active state, giving rise to MSW 
resonances (both in the neutrinos and anti-neutrinos channels). 

In the region of Am^ < 10^^ eV^, the production starts too late, namely after neutrino freeze-out. The only effect 
of the oscillations is then to redistribute the energy density between the active and sterile species, keeping the total 
entropy constant (there can be no "refill" from the thermal bath) so the bound on N^^^ does not apply. 

Let us then move to discuss the BBN probes ("'He and D abundances): at large Am^ their isolines essentially 
coincide with N^^^ . On the contrary, in the region of /S.m? < 10^^ eV^, if Ve — *■ oscillations occur then i/s 
are created by depleting u^, as just discussed. Thus the n/p ratio is affected and consequently the ^He abundance 
and, to a lesser extent, the D abundance. This is apparent in the i^i / Vs and V2 / Vs plots of fig. |3] where the bound 
extends in the lower part. The bound is stronger for the eigenstate which contains the larger portion of electron 
flavor: vi. In the case of v^/vs mixing nothing happens because no significant Ve component is present in v^, as 
013 is very small. 

In turn, for Am^ < 10"^ eV^ oscillations begin really too late, namely even after the decoupling of n ^ p re- 
actions. At that stage, the relative n/p abundance is no more affected by the neutrino populations (but only by 
neutron decay), so no bounds apply. 

The mixing with flavor eigenstates z^e, t'/j, Vr (lower row of Fig.|3} is qualitatively different, for the general reasons 
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Figure 4: The LSND anomaly interpreted as oscillations of 3+1 neutrinos versus the cosmological constraints. 
Shaded region: suggested at 99% C.L. by LSND. Black dotted line: 99% C.L. global constraint from other neutrino 
experiments (mainly Karmen, Bugey, SK, CDHS). Continuos red line: the bound from BBN corresponding to 
Yp ~ 0.26, i.e. = S.8 thennalized neutrinos. Dot-dashed orange line: the bound from LSS corresponding to 
nji^ = 0.01. 

explained in section|2] Essentially, since the flavor eigenstate is spread on two or three mass eigenstates, the active 
sterile oscillations occur at two or three different Arn^, of which one is always large enough to give an effect 
for every Amj4 (our variable of choice for the vertical axis). 

In the case of Ve/vs mixing, we see the effect of the solar mass splitting as a bump in the corresponding panel in 
fig-E] if the (mainly) sterile state lies below V2, the Arn^^j^ ^ 7 10~^ eV^ becomes the dominant mass difference.^) 
The bounds then approximate to vertical lines and they are stronger in tan 6s because of the resonant disposition. 
More precisely, this is true for the moderate effect expressed by the = 3.2 line; Am^^jj is sufficient to cause a 
more incisive effect {N^ = 3.8 or more) only in the case of ^He. 

In the case of v^^r/vg mixing, the Am^jj^^ plays exactly the same role, and this time it is large enough to affect all 
probes. 

Finally, the LSS structure bound on VLyh? consists of an horizontal line in the central part of the plots: qualitatively, 
as long as the sterile species is fully populated it is its contribution to the sum of the masses (i.e. the /S.m?) which 
is bounded from above. The constraints get weaker for very small tan 9s, where the sterile neutrinos are less 
efficiently produced {p^^ <C). At tan 9s > 1 the bound from fti, holds even for very small mixing, 9s ~ 7r/2, just 
because this region corresponds to heavy active neutrinos. 

In summary, fig.|3]displays the excluded, the allowed and the future testable regions of the active/sterile mixing 
parameter space, for six limiting cases, for what concerns cosmological probes. Any specific model of a sterile 
neutrino identifies a preferred point (or area) in one of these spaces or in a suitable combination of them, and 
should therefore be tested on them. This is what we do, as an example, for the LSND sterile neutrino in the next 
Section. 




3.2 LSND: in or out? 

The LSND experiment reports a signal |6| for i/^ oscillations in the appearance of De in an originally Df^ 

beam. The best fit point is located at sin^ 6'lsnd = 3 • 10"^, AmLgjMD = 1.2 eV^ (the whole allowed region is 
represented in fig.E). 

The "3+1 sterile" neutrino explanation assumes that the D^^ — > oscillation proceeds through Ufj^ ^ ^ i^e- 
The large LSND mass scale Atolsnd separates the three active states (split by the solar and atmospheric gaps) 
from the additional sterile state. The effective angle of the LSND oscillation 6'lsnd can be expressed in terms of 

^^Keep in mind that we are assuming a negligible 613 so there is no electron component in the third mass eigenstates which can feel the 
mixing with the sterile neutrino. That's why Am^^^ has no role in this case. 
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the two active-sterile angles 6 eg, dpis as 6'lsnd ~ 6es ■ dpLs-^^ However, each one of these two angles is constrained 
by several other experiments that found no evidence of electron or muon neutrino disappearance. Moreover, the 
— > Ve oscillations are directly excluded by KARMEN in a fraction of the parameter space. As a result 1221 . 
a large portion, but not the totality, of the area indicated by the LSND experiment is ruled out (see fig. 1^. The 
poor compatibility with solar and atmospheric oscillation data, in the context of four v mixing, also puts the LSND 
sterile neutrino in a difficult position 1231 . 

How does the LSND signal compare to the cosmological bounds discussed above? Fig. |3 shows the constraints 
from BBN and from LSS superimposed to the LSND region We see that the entire LSND region is ruled out 
by the BBN constraint: basically, for every value of its mixing parameters the LSND sterile neutrino completely 
thermalizes and implies an unacceptable modification of the *He primordial abundance. Of course, remember 
that allowing non-standard modifications to cosmology/BBN has the power to relax the bound to some extent. 
Exemplar is the case of a large primordial neutrino asymmetry, as discussed above. Other recent suggestions 
include l24l . 

The constraint on Vlp is well approximated by the horizontal line corresponding to Vl^h^ = m4/93.5 eV. It starts 
to bend (as discussed in Sec B.U only at smaller values of the effective 6'lsnd mixing angle. 

Fig. 15 that combines the precise computation of the evolution of mixed neutrinos in the Early Universe and of the 
neutrino experimental data, therefore reproduces and completes the estimates already presented in 1251 and 1261 . 

4 Sterile Neutrinos in Supernova sauce 

Generalities: Supernovse can be powerful and important laboratories for neutrino physics, and therefore in partic- 
ular for the physics of sterile neutrinos. The fundamental reasons for this basic fact are simply listed: 

• SNe are abundant sources of neutrinos, since this is the main channel into which most of their enormous 
energy is emitted; as a consequence, neutrinos play a crucial role in the evolution of the SN phenomenon; 

• given the characteristic temperatures of the SN environment, the typical energy of the emitted neutrinos 
(~ 10 20 MeV) is such that they can be easily detected on Earth; 

• SNe are so far away that neutrinos must travel over distances so large that they have plenty of time (or space) 
to experience and fully develop the consequences of several "exotic" effects (oscillations, non-conventional 
very feable interactions, decay...), if any is present; 

• SN cores are so extremely dense that neutrinos remain trapped and undergo matter effects that cannot be 
relevant anywhere else. 

On the other hand, it is true that the physics of supernovas is very complicated and demanding, and could pose a 
threat on their possible usefulness as "clean experiments". Nevertheless, the basic features are robust enough to 
be used as incontrovertible criteria, sometimes, maybe, requiring a sensible compromise between detailness and 
usefulness in the treatment of SN physics. At the present epoch, the twenty-something events of the SN1987a 
signal B7| already allow to set cautious constraints. The future is however brighter: running solar neutrino ex- 
periments could detect thousands of events from a future SN exploding at distance ~ 10 kpc and an even 
more impressive harvest of data could come from a future Mton water-Cerenkov detector or from other more SN- 
oriented future projects II28I . making the quantitative analysis of the SN neutrinos worthwhile for the search of 
sterile states. Previous analysis have investigated the subject II29I . 

Neutrino evolution: What we need to do is simply said: we must follow the fate of the neutrinos emitted from 
neutrino-spheres^' along their travel through the star mantle, the vacuum and (possibly) the Earth. ^' The existence 
of the sterile neutrino can introduce modifications at each of these steps, via matter or vacuum conversions, in 
different ways for each possible choice of the sterile mixing parameters. 

^^In our parameterization, this simply corresponds to a unit vector n = (ne,n^,nT) ~ 

'^'in constructing fig.|4]the BBN constraint lias been minimized setting Gea ~ f/js ~ ^^lsnd. when this is allowed by the neutrino 
disappearance data. 

Neutrino-sphere: the region of the star mantle at which the density becomes low enough that neutrinos (produced in the core) are no more 
trapped and freely stream outwards. 

'''For the range of Am^ < few eV to which we are confined by the cosmological bounds discussed in Sec|3] the MSW resonances with 
the sterile state occur out of the neutrino-sphere. The resonances would enter in the neutrino-spheres (in the inner core) for Am^ 10''' eV^ 
( 10*^ eV^ respectively). 
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In more detail: one has to follow the evolution of the 4x4 neutrino density matrix grn, written in the basis of instantaneous mass eigenstates. 
For instance, a Ue with energy E,y is described by Qrn, = VX, ■ diag (1, 0, 0, 0) ■ Vm where Vm depends on Ei, and, in general, on the position 
in the star. The mixing matrices in matter (Vm) and vacuum (V) are computed diagonalizing the Hamiltonian 

H= ^ [ydiag{m?,m|,m§,m|)yt +S^diag(Ve,V^,VV,0)l (16) 

and ordering the eigenstates according to their eigenvalues Hi = . /2E,y. Umi (Mm.4) is the lightest (heaviest) neutrino mass eigenstate 
in matter. The evolution up to the detection point is described by a 4 X 4 unitary evolution matrix U so that at detection point the density matrix 
g in the basis of flavor eigenstates is 

g = V -U ■ QmiEv) -U^ ■ with W = WEarth ■ ^vacuum ■ Wstar- (17) 

The evolution in vacuum is simply given by Wvacuum = diag cxp(~iLTn'^. /2E,j). Combined with average over neutrino energy it suppresses 
the off-diagonal element when the phase differences among eigenstates i and j are large. 

The evolution in the matter of the star is more complicated because several level crossings can occur, say at radii ri . . . rjv- At each one of 
them, there is a certain "jump" probability. In the present formalism, this can be expressed by a 4 X 4 rotation matrix P, with the rotation angle 
given by tan^ a = Pq /(I — Pc) where Pq is the level crossing probability. Indeed, in particular, when levels i and j cross in an adiabatic 
way, P = n. If instead the level crossing is fully non adiabatic P is a rotation with angle a = 90° in the (ij) plane. 

The computation of Pq at the crossings requires attention. Focussing e.g. on 6s < 7r/4, at a crossing between a mainly active state Ua and the 
mainly sterile state Us we compute it as 

Pc = = 7 = ^ = 7 ; ^ where sin 6*™ = n X sin ■ (18) 

e-y-1 ' dHa/dr ' 27r|cos26»™| as a 

where it is important to notice that 7 and 0™ must be computed around the resonance, where Tiaa = Ti-ss (or around the point where 
adiabaticity is maximally violated, in cases where there is no resonance) and are in general different from their vacuum values (that are instead 
conveniently used to parameterize Pq in the simpler 2v case). 1301 

Between a level crossing and the following one the evolution proceeds as governed by the matter Hamiltonian, so that the complete form for 

^^star is 

f2 m?, .\ I r^i 



Wstar = Pr„ ■ ■ ■ Pra ' diag cxp \ j ~2E^ j ^''^ '''^^ '^^^ } '^^ 2E^ ' ' 

In practice, given the importance of the matter effects in the star mantle and the very long baseline to Earth, the evolution between the level 
crossings and in vacuum averages to zero the off-diagonal elements that are possibly produced by the rotation matrices, introducing significant 
simplifications. However, if two states have Am^ < 10^^* eV^ vacuum oscillations do not give large phases: evolution in the outer region of 
the SN and in vacuum must be described keeping the off-diagonal components of the density matrix. 

Finally, for simplicity we now assume that the neutrinos do not travel through the Earth matter, so that WEarth is trivial. In the case of the 
SN 1987a signal, this was not true but disregarding it only implies a few percent eiTor, well within the general uncertainties for our purposes. In 
the case of the next SN event, it could be in general reintroduced. 

Having described the general formalism, let us now focus more closely on the peculiarities of the SN case. At the neutrinosphere, matter effects 
are dominant, so that matter eigenstates coincide with flavor eigenstates (up to a trivial permutation dictated by the MSW potentials that will 
be more clear below): the initial density matrix consists of Qm = diag {^%^ , 'I'^^ , , ^1^)/^tot' where are the flavor fluxes from the 
neutrinosphere and ^tot stands for their sum. The final fluxes at the detection point will be then given by the diagonal entries of g of eq. 1171 : 
= diag(e) ■ $tot. 

Concerning the initial fluxes, the accurate results of simulations are usually empirically approximated by a Fermi-Dirac spectrum for each 
flavor Ue,Ue, Vx {Vx collectively denotes i^fi,T, ^ii,t), with a "pinching" that slightly suppresses the low energy portion and the high energy 
tail. Based on the recent results of 1311 , we adopt the following average energies and total luminosities for the various neutrino components 
at the time of the snapshot of fig. |^ (see below): {E^j^^v^ ) ~ 12, 14, 14MeV, L„^, v^,u^ — 30, 30, 20 ■ 10"* "'^ erg sec ^, assuming (in 
accordance with numerical calculations) that the ratios of luminosities do not vary much during the whole emission. The initial flux of sterile 
neutrinos is assumed to be vanishing, as a consequence of the fact that matter oscillations only take place out of the neutrinosphere. 

The MSW potentials of eq. 1161 , experienced by the neutrinos in SN matter are 

Ve = \/2Gpns (SYe - 1) /2, V-r = V^ + V^^, 

= V2GFnB (Y, ~ 1) /2, 14 = 0, 

where is the baryon number density (rig = p{r) /mjv where mjy ~ 939 MeV is the nucleon mass ) and Ye = {N^- — )/ng is the 
electron fraction per baryon. Antineutrinos experience the same potentials with opposite sign. The difference V^t in the and Ur potentials, 
which appears at one loop level due to the different masses of the muon and tau leptons 1331 . is, according to the SM 
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The effect is not irrelevant in the inner dense regions: for densities above p ~ 10** gem""', the fir vacuum mixing is suppressed. 

A crucial point concerns the characteristic of the matter density p(r) and of the electron fraction Ye in the mantle of the star. We adopt the 
profiles represented in fig. Island we model them with analytic functions that preserve their main features. Namely, the density profile decreases 
according to a power law r~'^ out of the ~ 10 km inner core (where instead it has a roughly constant, nuclear density value). At much larger 
distances the density profile gets modified in a time-dependent way by the passage of the shock wave. Present simulations have difficulties 
in reproducing this phenomenon and therefore cannot reliably predict the density profile in the outer region. Therefore for r ^ 500 km we 
assume a power law p(r) = 1.5 1O'*(P0 /r)'^g/cm^, which roughly describes the static progenitor star. 

In turn, the peculiar Ye profile in fig.|5]is inevitably dictated by the deleptonization process: behind the shock wave which has passed in the 
mantle matter, the electron capture on the newly liberated protons is rapid, driving Ye to low values (~ 1/ 4). In the outer region, where 
the density is sensibly lower, the efficiency of the capture is much lower, so that Ye essentially maintains the value ~ 1/2 typical of normal 
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Figure 5: Inside the supernova. Left panel: Density p{r) and electron fraction Yf^{r) from H32V . Center panel: 
matter potentials for ( red solid line) and for v^^t (blue solid line); the dashed lines are the analytic modelization 
that we adopt. Right panel: a representation of the antineutrino matter eigenstates with their flavor content, as 
functions of the radius in the interior of the SN; the drawing is made in the specific case of small Ve/vs mixing and 
large I\m\^{> I^m^^^) for the sake of illustration; we are assuming small Oi^ and normal hierarchy. 

matter. ^''^ This is important because the matter potential Ve of electron (anti)neutrinos flips sign, see eq. 1201 , when, in the deep region of the 
mantle Ye steeply decreases below 1/3. 

The knowledge of the matter potentials (and of the neutrino masses and mixings) allows to draw the pattern of the (anti-)neutrino eigenstates 
with their flavor content, as functions of the radius r in the interior of the star. This is represented in the right panel of fig.|5] in the specific case 
of small Ue/vs mixing and large Am^(2> Am^^j^j) for the sake of illustration. 

In the end, we collect the modified (anti)neutrino fluxes and deduce the modifications to some relevant observables; 
in particular, we focus on the final flux of v^, which are best detected through v^p e+n at the Cerenkov detectors. 



4.1 Results 

The results are collected in fig.|6| where we plot the reduction of the i>e event rate in a typical Cerenkov detector 
due to sterile mixing and in fig-El where we plot the modified average energy and the distortion of the Dp 
spectrum for a few specific choices of mixing parameters. 

In order to understand qualitatively the main features of these results it is useful to look at the pattern of level 
crossings, like the one qualitatively depicted in fig. El Indeed, as discussed above, the active/sterile MSW reso- 
nances in the matter of the star are the crucial places where the neutrino flux is non-trivially modified, the rest is 
simply vacuum oscillations. There are three possible resonances: 

1. The mostly Dg eigenstate crosses the mostly i>e eigenstate at r ^ 100 km, where Ve flips sign. At this 
point matter effects dominate over active neutrino masses, so that active mass eigenstates coincide with 
flavor eigenstates. Since Ve flips sign in a steep way this resonance is effective only if Amf^ > 10^ ^^'^ eV^ 
(different SN simulations gives values in this range). 

2. If the mostly sterile eigenstate is the lightest one (in our parameterization this needs 6's > 7r/4) the two 
eigenstates in 1 . cross again at larger r. Pictorially, this second resonance is present when the sterile black 
line is lowered below the others in fig.|5] This MSW resonance occurs at large r where Ve is smooth, so that 
it is effective down to Am^^ > 10~^"^^ eV^. Again, the significant uncertainty is due to uncertainties on the 
SN density gradient. 

3. When the mostly sterile eigenstate is the heaviest or the next-to-heaviest state, it also crosses one or both of 
two mostly iy^^r eigenstates. This is the case illustrated in fig.|5] The values of Ani^^ and Am§4 determine 
at which r these crossings takes place, and consequently the flavor composition of the mostly active states at 
the resonance. Entering in the SN, the small Dg component of disappears as soon as 14 — dominates 
over ATOgmj. In any case, active/sterile MSW resonances with the mostly D^,r states affect only marginally 
the De rate, right because the 3 contribution to De is secondary. 

^"^The data refer to ~0.3 sec after bounce for a typical star of ~ 11 solar masses. The subsequent evolution is supposed to move the wave of 
the Ye profile slightly outwards, maintaining, however, its characteristic shape. 
^^'We focus on Uep scatterings with the cuts and efficiency of the Kamiokandell experiment. The cross section is taken from 1341 . 
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At each of the resonances, part of the neutrino flux can then convert into the sterile state and deplete the final 
active flavors. Combinations of more than one resonance can occur, depending on the mixing parameters. These 
considerations allow to understand the regions of i>e flux reduction in fig.|6l 





Figure 6: Sterile effects in supernovas. The iso-contours correspond to a 10, 20, 30, 40, 50, 60, 70 % deficit of the 
SN Vf, total rate due to oscillations into sterile neutrinos. The deficit is measured with respect to the rate in absence 
of active/sterile oscillations but of course in presence of active/active oscillations (which reduce the unrealistic no- 
oscillation-at-all rate by ^ 10%j. We shaded as disfavoured by SN 1987a data regions with a deficit larger than 
70%. While the qualitative pattern is robust, regions with MSW resonances can shift by one order of magnitude in 
Am^ using different SN density profiles. / Dg mixing ( not plotted) does not give significant effects. Fig. ^studies 
in detail the sample points here marked as A, B, C, and the regions 1, 2, 12, 3 are discussed in the text. 

Let us start from the case of vi j Vg mixing. Resonance 1 gives a sizable reduction in region 1 and resonance 2 gives 
a sizable reduction in region 2. Had we ignored solar mixing the maximal deficit would have been 100%, while 
in presence of solar oscillations the maximal effect is a 80% deficit. More precisely, in the interior of region 
1 one obtains — sin^ because resonances 1 and 3 are fully adiabatic. In the interior of region 2 one 

obtains = sin^ Qsvl\\'^%^ ^ because resonance 2 is fully adiabatic and resonance 1 irrelevant. Therefore, given 
the assumed initial fluxes, the rate gets reduced slightly more strongly in region 2 than in region 1 . 
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Figure 7: Sterile effects in supernovse. Average z^p energy in the vi/vs plane (iav? 9^ > 1 portion) and distortion 
of the Deflux at sample points A, B, C. 

In region 12 both resonances 1 and 2 are effective, and tend to compensate among each other: resonance 1 converts 
Ve into Vs and resonance 2 reconverts i>s into Ve- 
in region 3, resonances 3 gives a 20% suppression of the Ve rate, that sharply terminates when Amlj^ < Am^^^^. 
This is due to a strong suppression of the mostly t- eigenstates, which due to solar oscillations would give a 20% 
contribution to the z/g rate (ignoring solar mixing, there would be no suppression of the i>e rate in region 3). This 
reduction of the D^,,r fluxes induced by resonances 3 could be better probed by measuring the NC rate (which gets 
a < 40% reduction) and, if neutrinos cross the Earth, by distortions of the De energy spectrum. 
The tail at smaller Am^ and around maximal mixing is due to vacuum oscillations, that can reduce the rate by 
<50%: their effect persists down to Am^ _ e^^/D ~ 10"i**eV^ The precise value depends on the distance D 
(we assumed D = 10 kpc). 

The other mixing cases are understood in similar ways. We remind that our parametrization is discontinuous at 
9s = 7r/4: this is reflected in the panel of fig. |6l which illustrates 1^2/ I's mixing. Resonance 2 sharply terminates 
when Am24 < Amg^^. Except for these differences, this case is quite similar to the previous one, because i>i 
and D2 get strongly mixed by matter effects. On the contrary 1^3/ i^s mixing (not shown) does not give a significant 
reduction of the De rate. 

Mixing with the flavor eigenstates behaves in a similar way. Namely, the v^^^/vs figure shows the reduction 
in region 3 due to resonances 3, and the v^/vs figure shows the reduction in region 1 due to resonance 1. The 
additional feature at Am\^ ~ Ato^^jj is due to adiabatic conversion (for large sterile angles) between the mostly- 
sterile state and that are almost degenerate in this condition. For larger Amf 4, there is no component in 
so that the crossing is totally non adiabatic, while for smaller Am^^ the two states are separated. Sterile effects 
persists at all values of Am\j^, even if it is small. 
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Figure 8: Summary of sterile neutrino effects. The shaded region is excluded at 99% C.L. (2 dof) by solar or 
atmospheric or reactor or short base-line experiments. We shaded as excluded also regions where sterile neutrinos 
suppress the SN1987A rate by more than 70%. This rate is suppressed by more than 20% inside the dashed blue 
line, that can be explored at the next SN explosion if it will be possible to understand the collapse well enough. 
Within standard cosmology, the region above the red continuous line is disfavoured (maybe already excluded) by 
BBNandLSS. 

The reductions highlighted in fig. |6lhave to be compared to the measured flux to set constraints and identify 
anomalies. This will certainly be fruitful in the future. At present, however, the theoretical uncertainties on the 
SN evolution and fluxes and the smallness of the SN1987a data set only allow to put conservative constraints: we 
simply shaded as 'disfavoured' regions where sterile effects reduce the rate by more than 70%. 

In the same perspective, important observables in the future SN event will be the spectrum, where distortions 
could be induced by sterile oscillations, or more generally the average Vp energy. In absence of sterile oscillations 
we expect {E^J) « 15MeV with a quasi-thermal spectrum. In presence of sterile oscillations, these observables 
can be affected in selected regions of the parameter space, as illustrated by figure The average energy can 
increase up to ~ 18MeV and decrease down to w 12MeV along the sides of the MSW triangle. The 

distorted spectrum which correspond to the first case is plotted in fig. ^Jj. Vacuum oscillations can also give 
well known distortions, as exemplified in fig. Et- These effects seem larger than experimental and theoretical 
uncertainties. In all other cases sterile effects give a quasi-energy-independent suppression of the i>e rate, and 
therefore negligibly affect {EpJ. Fig.0\ gives an example of this situation. 

5 Conclusions 

We have performed in (Q a systematical study of the effect of an extra sterile neutrino in all possible contexts, for 
any choice of its mixing parameters with the active neutrinos and fully including the active/active mixings now 
established by the results on solar and atmospheric oscillations. We considered cosmology (BBN, CMB, LSS), 
astrophysics (the Sun, SN...) and terrestrial neutrino experiments (atmospheric neutrinos, reactors, short- and long- 
base line beams). In these Proceedings I presented in detail the analysis relative to the effects in cosmology and in 
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supernovse focussing on the principles of the computational techniques and on the understanding of the results. 

We find no evidence for a sterile neutrino in the present data; figure |8l collects the presents constraints, in the 
six limiting cases of mixing with a mass eigenstate j^i,2.3 or a flavor eigenstate i^e.^i.r- Every specific model of 
a sterile neutrino identifies a point or a region in one of these spaces or in a suitable combination of them, and 
should therefore be compared with the reported bounds. In particular, 1 showed in Sec l3.2l that the region of the 
LSND sterile neutrino is excluded by the cosmological constraints (unless they are relaxed by some non-standard 
modification to cosmology). 

Finally, 1 discussed promising future probes of sterile effects. In the context of the cosmological observations, it 
looks important to improve the measurements of the primordial ^He and Deuterium abundances, overcoming the 
systematic uncertainties. It looks also especially important to exploit the complementarity of the two probes, that 
are differently sensitive to sterile neutrinos in several regions of the parameter space, also in order to constrain the 
non-standard cosmological modifications. Moreover, future measurements of the CMB and LSS power spectra 
will expand the tested range of sterile parameters. Concerning supernovas the best improvement would come from 
. . . the occurring of a new explosion, which would allow to probe sterile effects through the predictions on the 
neutrino flux and spectrum, in regions that are not easily accessible to other tools (thanks to the long base-Une 
and the extreme matter effects). Although progress must also come on the overall theoretical uncertainties of SN 
models, most of the results could be neat enough. 

Other important probes for sterile neutrino effects are discussed in (!]: among them, solar neutrino experiments at 
sub-MeV energies and proposed reactor and LBL experiments look most interesting. 
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